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Abstract 

We investigate the thermodynamics of a crystalline solid applying q-deformed algebra of Fibonacci os¬ 
cillators through the generalized Fibonacci sequence of two real and independent deformation parameters 
qi and q^. We based part of our study on both Einstein and Debye models, exploring primarily (q\, (p)- 
deformed thermal and electric conductivities as a function of Debye specific heat. The results revealed that 
^-deformation acts as a factor of disorder or impurity, modifying the characteristics of a crystalline struc¬ 
ture. Specially, one may find the possibility of adjusting the Fibonacci oscillators to describe the change of 
thermal and electrical conductivities of a given element as one inserts impurities. Each parameter can be 
associated to different types of deformations such as disorders and impurities. 
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I. INTRODUCTION 


The interaction between atoms allows propagation of elastic waves in the solid medium which 
can be both transverse and longitudinal. If the oscillations of the atoms around the equilibrium 
positions are small, which should occur at lower temperatures, the potential energy of interaction 
can be approximated by a quadratic form of the displacements of atoms from their equilibrium 
positions. A crystalline solid, whose atoms interact according to this potential, is called a harmonic 
solid. In harmonic solids, elastic waves are harmonics and the normal modes of vibration in 
crystalline solids [IQ. A large number of phenomena involve quantum mechanical motion, in 


particular thermally-activated particles, obeying the T 3 law. Thermal excitations in the system are 
responsible for phonon excitation [2L13]. 

The study conducted by Anderson, Lee and Elliot [4-(|] shows that the presence of defects 
or impurities in a crystal modifies the electrostatic potential in their neighborhood, breaking the 
translational symmetry of the periodic potential. This perturbation can produce electronic wave 
functions located near the impurity, ceasing to be propagated throughout the crystal. 

A possible way to generate a deformed version of the classical statistical mechanics consists in 
replacing the Gibbs-Boltzmann distribution by a deformed version. In this respect it is postulated 


a form of deformed entropy [7] which implies a generalized theory of t 
We apply the g-deformation in the models of Einstein and Debye [8 


lermodynamics. 


ION , and our results show 


that the factor q acts as an impurity, modifying the thermodynamic quantities such as entropy, 
specific heat, thermal conductivity, etc. 


In t 


lis work we insert the given parameters of deformation q x and q 2 , called Fibonacci oscil- 


in [[13 


15 


lators [121, which is a formalism recently proposed in the (-/-calculation that has been investigated 


]. They provide a unification of quantum oscillators with quantum groups, keeping the 


degeneration property of the spectrum invariant under the symmetries of the quantum group [16- 


201 . The quantum algebra with two deformation parameters may have a greater flexibility when 


it comes to applications in realistic phenomenological physical models 112 lL 1221 and may increase 
interest in physical applications. 

The paper is organized as follows. In Sec. [Ill we present the deformed algebra. In Sec. [Till 
we apply the Fibonacci oscillators, in the Einstein IIII Al and Debye IIII Bl models, and finally, in 
Sec. tlTVl) we make our final comments. 
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II. ALGEBRA OF THE FIBONACCI OSCILLATORS 


It is well-known that the generalization of integers in general is given by a sequence. A basic 
procedure in g-algebra [11] is a generalization of integers. Two well-known ways to describe a se¬ 
quence are the arithmetic and geometric progressions. A simple generalization that encompasses 
both of is the Fibonacci sequence, which as we know is a linear combination where the third 
number is the sum of two predecessors, and so on. Here, the numbers are in that sequence of gen¬ 
eralized Fibonacci oscillators, where new parameters (q \, q 2 ) are introduced. Thus, the generalized 
spectrum may be given by the whole Fibonacci number. 

The algebraic symmetry of the quantum oscillator is defined by the Heisenberg algebra in terms 


of the annihilation and creation operators c, e ', respectively, and the number operator N [113L 
via 


230 


t 2 t 2 n; 

c * c i - <h c i c i = 42 


[N, c+] = c\ 


and 


J 


,2 t 


^271.7, 


[N,c] = -c. 


(1) 

( 2 ) 


In addition, the operators obey the relations 


Jc=[N], <x? = [l + N], 


,t _ 


(3) 


[1 + ni,q 1 ,q 2 ] — 4l \ n i,qi,qi\ + 42 % ° r [1 3“ n i,qi,qz\ ~ 42 [^*,91, 92 ] 3~ C(i 


2n,; 


(4) 


The oscillator [16 


171] allows us to write the (q u g 2 )-deformed Hamiltonian [2Q0 as follows 


U = c,J\. 


(5) 


The Fibonacci basic number is defined by [24] 


[ n i,qi,q2] C i C i ~~ 


2 tii 2 rii 

42 - 4i 

9 9 ' 

42 - 4i 


( 6 ) 


where q\ and q 2 are real positive and independent parameters of deformation. 
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III. APPLICATION OF FIBONACCI OSCILLATORS 


A. (q\ , g 2 )-deformed Einstein solid 


We consider the solid in contact with a thermal reservoir at temperature T, where we have re¬ 
labeling the j-th oscillator Given a microscopic state { n 3 } = {n \ , n 2 ...., un}, the energy of this 
state can be written as, 

E i n j} = 

where u>e is the Einstein frequency characteristic. We can obtain ( r/i. g 2 )-deformed energies from 
the definition of the Hamiltonian ©, and the definitions provided earlier, 


E, n , 


hhJEfr i . r . . hw E (2 ln(g 2 ) - 2 ln(gi)) n 

= — (fa, 91 , 92 } + fa, 71,72 + 1 ]) = — + -^ 3^2 -’ ( g ) 


and when qi = q 2 = 1, we recover the usual spectrum 


E n = (2 n + 1). 


(9) 


With the result of the Eq.®, we can rewrite the partition function in the form, 


N 


"( 71 , 72 ) 


^ exp [~PE n . qi J 


7 n 

"( 1 , 71 , 92 )’ 


n=0 


( 10 ) 


where 


1 ( 1 ,<7l ,<?2) 


ql a 


exp ( a(q \ 2 ) (q\ a — q% a ) 

We define a (q 1: q 2 )-deformed Einstein function E(a) qijq2 , 


(ID 


E(Ot) q i,g 2 


f 2aq 1 q 2 [ln(qi) - ln(g 2 )] 

V fa 2 " - qt) 


where a = 


©/ 


T(ql - ql) 


( 12 ) 


As one knows ©^ is the Einstein temperature, defined by 


©K = 


HuJe 

Kb 


(13) 
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and kb is the Boltzmann constant. When q\ = q 2 —>■ 1, we have the undeformed function 


E(a*) 


exp (a*)a 2 * 
(exp a* — l) 2 ’ 


where 



( 14 ) 


We can determine (qi, r/ 2 )-dcformed the Helmholtz free energy per oscillator and entropy, re¬ 
spectively 


/, 


91 >92 



lim — In 5 

N—too N 


( 91 , 92 ) 


—KbT In 


exp 


q\ a 


»(9?-9? 2 ) 

2 


(qf a 



(15) 


S', 


91,92 


d/91,92 

dT 


= k b < In 


q\ a 


exp 


2aql a (\w(q 1 ) - ln(g 2 )) 

0 qT ~ ql a ) 


+ 


®(ql - qi) 


ql a )\ 


+ 


(16) 


In Fig. ([[]) is shown the behavior of the entropy as a function of temperature variation. We 
observe that all the curves have the same behavior at low temperature. However, as the temperature 
increases the role of the (/-deformation becomes much more evident. For instance, note that q 2 
tends do decrease the entropy more than q x . As we anticipated, these parameters can play different 
roles. While one can affect disorders the other may control impurities. 

Now we determine (q 1 , g 2 )-deformed specific heat, and we can do it by inserting the Einstein 
function E(a) qitq2 , defined by Eq.(fT2l). into equation below 


<v„.jr) = r 


as, 


91,92 


dT 


— kb 


f 2a<?i(/ 2 [ln((/i) - ln(g 2 )] 

V (■ ql a " ' 


ql a ) 


(17) 


%«( t ) = 3k ^H 1,92- (18) 

The complete behavior is depicted in Fig. ©. One should note that when T (~) E , and thus the 
ratio a = y < 1, and @g around 100A" for common crystals, one recovers the classical result 
c v q 1 n —» 3 kb, known as the Dulong-Petit law. However, for sufficiently low temperatures, where 
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- ql=0.1, q2->0 ql=0.5, q2^0 - ql=l, q2^0 

-ql=0.1, q2—»1 - ql=0.5, q2->l — ql=l, q2^1 


- ql=0.1, q2—»0 ql=0.5, q2->0 - ql=l, q2->0 

-ql=0.1, q2—>1 - q 1=0.5, q2->l ql=l, q2^1 



Figure 1: ^i, ^-deformed entropy S qii q 2 vs temperature T in the following intervals: (T=0, ••• , 10 K) (left) 
and (T=0, • • • , 1000FQ (right) 


- ql=0.1, q2->0 ql=0.5,q2->0 ■ ql=l,q2->0 

-ql=0.1g2—>1 - ql=0.5,q2—*1 ql=l,q2^1 




Temperature T (K) 


- ql=0 

1, q2—0 

ql=0.5, q2-0 - 

ql=l 

,q2^0 

— ■ •ql=0 

1, q2—4 

- ql=0.5, q2—*1 

ql=l 

,q2-4 



Figure 2: q\, ^-deformed specific heat cy n , l2 vs temperature T in the following intervals: (T=0, • • • , 1 A') 
(left), (T=0, • • • , 10A') (center) and (T=(), • • • , 1000 A') (right). 


T -C Qe and therefore a»l, specific heat decreases exponentially with temperature [2], as 


Cy —>■ Kb 




(19) 


In general, the invariance of specific heat at high temperatures and its decrease at low tem¬ 
peratures show that the Einstein model is in agreement with experimental results. However, at 
sufficiently low temperatures, specific heat does not experimentally follow the exponential func¬ 
tion given in Eq.(fl9l). As for the (V/i. (^-deformed case we see a significant change in the curves 
at intermediate temperatures. 
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B. (q\, g 2 )-deformed Debye solid 


Corrections of Einstein model are given by the Debye model, allowing us to integrate from 
a continuous spectrum of frequencies up to the Debye frequency tj D , giving the total number of 
normal modes of vibration [1-3] 


g(u)du = 3iV, (20) 

where g(u)du denotes the number of normal modes of vibration whose frequency is in the range 
(c u, u + du). The function g{u), can be given in terms of the Rayleigh expression as follows 



8tt 




u 2 duj 
7 r 2 c 3 ’ 


( 21 ) 


where c is the speed of light and A wavelength. The expected energy value of the Planck oscillator 
with frequency u s is 


(Es) 


hu s 



( 22 ) 


Using Eqs.(l2TI) and (l22l ) . we obtain the energy density associated with the frequency range (u, cu + 

du). 


u(u)du 


h u 3 du 

(*)-f 


(23) 


To obtain the number of photons between u and u + du, one makes use of the volume of the region 
on the phase space [2], which results in 


g(u)du ~ 


2V 

1J 



Vu 2 du 
7 r 2 c 3 


(24) 


Thus, replacing Eq.([24l) into Eq.([2TT). we can write the specific heat for any temperature. We now 
apply (qi, g 2 )-deformation in the same way as in Eq.(fT8l). 


°v qi , q2 ( T ) - 3«B-D(ao, llW ), 


(25) 
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where D{olq ) is the (gi, g 2 )-deformed Debye function, defined by 


D{ a 0q 1: q 2 ) 


" a °‘H,<i2 Q^exofct) huin On 

LX CAp _ Uqi,q 2 1,12 / O /T\ 

UjOL) ^0 ai no m rri 1 


K 1 ,J 3 ./o [exp (ct) — 1] 


KrT 


T 


ud 0 i - 2wd 


ln(g 2 ) - ln(gi) 
?2 - 0.1 


(27) 


where cu£> and 0 Dq q , are the (gi ■ g 2 )-deformed Debye frequency and temperature, respec¬ 
tively, and Wfl is the Debye frequency characteristic. Integrating Eq.(l26l) by parts one finds 


■^( Q;0 91, 92 ) 


Scrn 


ex P(«0 9l , 82 ) - 1 


+ 


12 


f«0, 


f, i -"2 a 3 da 


a; 


o 91 , 92 7o exp (a) — 1 


(28) 


which can be integrated out to give the full expression 


, 99 ) 


3a ( 




+ 


12 


exp(a 0 ) — 1 ' «n 1 15 4 

V U 91 >92 2 Uqi , Q 1 


7T 4 & 0 o 

-_)_ a 3 ln[l — exp(a 0 )] 

K /I U 91,<22 1 U 91,92'J 


+ ^On.aL^fexplao,!^)] - 6ao 91ira Li 3 [exp(ao, ll , 2 )] + 6Li 4 [exp(a 091 , 92 )] >, (29) 


where 


U n{z) = Y J hi 


k =0 


(30) 


is the polylogarithm function. For T 0 Dq , a 0 ( gi , g2 ) ^ 1, then the function D(a 0 ) q can be 
expressed in a power series in a 0 


(Xn 


D(a 0 ) q - 1 - + 


(31) 


so that for 


T —> OO , CVq —y 3 txB- 


(32) 


On the other hand, for T -C 0n q , «o 9 1, then we can write function D(a 0 ) q as 


a, 


12 a 3 da ^ 

— / -^ + O[exp(-a 0 J], 


exp (a) — 1 


(33) 
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47r 4 

5 «L 


47T 4 


T 


( 34 ) 


Thus, as in the usual Debye solid, the low-temperature specific heat in a q-deformed Debye 
solid is proportional to T 3 , due to phonon excitation, a fact that is in agreement with experiments. 
Thus, let us express the (< 71 , g 2 )-deformed specific heat for low temperatures as follows: 


,92 


127T 4 7tB 


T 


= 1944 


Of) 1 \0 D 


T 


J 


molK 


(35) 


For the (q 1 , g 2 )-deformed case one can observe the changes that occur with Debye temperature, 


specific heat, thermal and e 


conductivity (ac) — see [25], we obtain 


ectrical conductivies. By using the relationship established for thermal 


« = - C V vl , 


(36) 


where v is the average velocity of the particle, Cy is the molar heat capacity and l is the space 
between particles. We can deduce a relationship between the thermal and electrical (a) conduc¬ 
tivities through the elimination of l (as a = where m is the electron mass, n is the number of 
electrons per volume unit and e is the electron charge), such that 


K 1 Cymv' 2 

a 3 ne 2 


(37) 


In a classical gas the average energy of a particle is | mv 2 = § k b T, whereas the heat capacity is 


so that 


K 


= U»lYt. 

a 2 V e / 


(38) 


The ratio is called the Lorenz number and should be a constant, independent of the temperature 
and the scattering mechanism. This is the famous Wiedemann-Franz law, which is often well sat¬ 
isfied experimentally, and the Lorenz number correctly given [25]. By using the [q \, <y 2 )-deformed 
relations presented above, we start from Eqs. (l36l) and (1381) to determine the important relations for 
(< 71 , q 2 )-deformed thermal and electrical conductivities 


K qi, <72 — 


KC Va 


91.92 


C V 


and 


a q{ -<j 2 “ 


K 


(39) 
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Recall that to compute these deformed quantities in terms of the specific heat cy q q we make 
use of Eq.([29l) and its suitable limits. We present in Tab. (IB, changes that occur with the De¬ 
bye temperature, specific heat, thermal conductivity and electrical conductivity, of some chemical 
elements. 


Element 

U Q 

7^) 

C v q 

Kq 

a™ 

Oq 


<?2=1 

92=0.5 

92=0.1 

92=1 

92=0.5 

92=0.1 

92=1 

92=0.5 

92=0.1 

92=1 

92=0.5 

-Q 

to 

II 

O 

Pb 

105 

194 

488 

4.53 xlO 4 

7.2x10^ 

450 

0.35 

0.06 

0.003 

0.48 

0.076 

0.00477 

Bi 

119 

220 

554 

3.11xl0 4 

4.9xl0 d 

309 

0.08 

0.013 

0.0008 

0.09 

0.0143 

0.00089 

Yb 

120 

222 

558 

3.04xl0 4 

4.8x10^ 

302 

0.35 

0.055 

0.0035 

0.38 

0.06 

0.0038 

Pt 

240 

444 

1116 

TsxitF - 

601 

38 

0.72 

0.11 

0.007 

0.96 

0.15 

0.01 

Pd 

274 

506 

1275 

2.55 xlO 15 

404 

25 

0.72 

0.11 

0.007 

0.95 

0.15 

0.01 

Y 

280 

518 

1302 

2.39X10 4 

379 

24 

0.17 

0.027 

0.002 

0.17 

0.027 

0.002 

Zn 

327 

604 

1521 

TJxicF 

238 

15 

1.16 

0.18 

0.01 

1.69 

0.27 

0.017 

Mn 

410 

758 

1907 

762 

121 

7.6 

0.08 

0.013 

0.0008 

0.072 

0.011 

0.0007 

Ti 

420 

776 

1954 

708 

112 

7 

0.46 

0.073 

0.005 

0.23 

0.04 

0.002 

Ni 

450 

832 

2093 

576 

91 

5.7 

0.91 

0.14 

0.009 

1.43 

0.23 

0.014 

Fe 

470 

869 

2186 

506 

80 

5 

0.80 

0.13 

0.008 

1.02 

0.016 

0.01 

Os 

500 

924 

2326 

420 

66 

4.2 

0.88 

0.14 

0.009 

1.10 

0.17 

0.011 

Ru 

600 

1109 

2791 

243 

38 

2.4 

1.17 

0.18 

0.01 

1.35 

0.22 

0.013 

Cr 

630 

1165 

2931 

210 

33 

2 

0.94 

0.15 

0.009 

0.78 

0.12 

0.0077 

Si 

645 

1192 

3000 

196 

31 

1.9 

1.48 

0.23 

0.015 

- 

- 

- 

Be 

1440 

2662 

6698 

18 

2.8 

0.17 

2.00 

0.32 

0.02 

3.08 

0.49 

0.03 

C 

2230 

4122 

10374 

4.7 

0.75 

0.048 

1.29 

0.2 

0.01 

- 

- 

- 


Table I: Chemical elements and their respective Debye temperatures “(if), Specific heat b ( —4^ ), thermal conduc¬ 
tivity c ^^, Electrical conductivity d (ohm ■ cm ) -1 x 10 5 , for T = 300if and the g-deformed values for q\ —>• 1 
and q -2 = 0.1, q 2 = 0.5 and <72 = 1. Ill] 


For illustration purposes, we choose iron (Fe) and chromium (Cr), two materials that can be 
employed in many areas of interest. In Figs.® |4]) we present deformed values of Fe ( Fe q -d e f ) 
(black) and Cr ( Cr q _ de f ) (green), for values q\ — 1 and q 2 — 0.1, ■ • ■ ,1, where for this range we 
assume the maximum deformation (q\ = 1 and q 2 = 0.1) and the pure element (bulk) (q\ = 1 and 
q 2 = 1). The other elements are represented by colors and indicated in the very figure. 

On the left side of the Fig.©, we can observe that before reaching their limits, black and green 
curves can assume the values of Debye temperatures (0 D ) of other elements. The Fe q _ de j e.g., 
equates to: beryllium (Be) when q 2 ~ 0.23, chromium (bulk) (Cr) q 2 & 0.75 and osmium (Os) 
q 2 r; 0.94. On the right, we have the behavior of the curves obtained for the specific heat cy . We 
note that the behavior is quite different from the previous curves d D , i.e., the curves start at lower 
values (maximum deformation) until they reach their pure values. Having Fe q -d e f as an example 
again, it is possible to see, as it reaches the value of specific heat capacity of all the elements, 
including Cr (bulk) when q 2 rj 0.74. 

On the left side of the Fig.©, we have the Cr q -d e f taking over values: manganese (Mn) when 
q 2 w 0.4, titanium (Ti) with q 2 « 0.77 and Fe bulk q 2 rs 0.95. On the right, we have the behavior 
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Figure 3: Debye Temperature Op depending on the variation q 2 = 0.1, • • • , 1 and q\ = 1, (left). Specific 
Heat cy depending on the variation q 2 = 0.1, • • • ,1 and qi = 1, (right) 


-Feq-def - Crq-def - Mn - Y Ti 

- Pt -Fe - Ni 


- Crq-def-Feq-def - Bi - Ti Yb 

- Pb - Cr - Pd 




Figure 4: Thermal Conductivity k depending on the variation q 2 = 0.1, • • • , 1 and q\ = 1, (left). Electrical 
Conductivity cr depending on the variation q 2 = 0.1, • • • ,1 and q\ = 1, (right) 


of the curves obtained for the electrical conductivity (cr), where we observe that the ytterbium (Yb) 
e.g., has its value reached by Fe q _ de f for q 2 ~ 0.7 and the Cr q _ de f q 2 ~ 0.78. Notice that in the 
present case q 2 develops an effect of impurity of the material. Such that the more q approaches 
zero the less thermal and electric conductivities approaches zero too. This is in accord with the 
experimental measures of conductivity of some good conductor that reduces e.g., its electrical 
conductivity by doping it with impurities. 

Let us now return to Eq. (l29l) . where we have the complete Debye function D(a 0 ) to see the 

behavior of the thermal conductivity. Thus, in the Fig.©, we show a comparison to the thermal 
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- q!=0.2,q2=l.8 - ql=l,q2=l - q1=1,q2=1.4| 



Temperature (K) 

Figure 5: Thermal conductivity k of Fe as a function of temperature in the range of T = 0, • • • , 1200 K, 
and values for q\ and q 2 given explicitly. 


conductivity n for a pure and impure material. We have the thermal conductivity as a function of 
temperature (T) for Fe (bulk) and a combination of the values q\ and q 2 for the Fe q -d e f (impure). 
Therefore, we have the Fe (red) when q 1 — q 2 — 1 (pure), and when q x = 0.2 and q 2 = 1.8, 
we have a similar behavior to silicon (Si) (golden), whose value for the Debye temperature is 
0d = 645A'. Finally, for a combination of values q\ — 1 and q 2 = 1.4 we have a curve that is 
similar to that of zinc (Zn) (orange), with d D = 327 K. 

We note that there are a large number of combinations of q\ and q 2 parameter values to be 
tested. In Tab. (|H) we show only two options (q\ = l,q 2 = 0.1) and (q 1 = 1 ,q 2 = 0.5) and 
Figs. (l3l |4]) it is shown that one pure material gets impurities by doping, for instance, it may present 
properties of other [IQ. The results of this study with two deforming parameters differ from the 
results previously obtained in [8fl by considering only one parameter. Now is clear there exists 
another one parameter that can play a different role of the other. 


IV. CONCLUSIONS 

We apply the Fibonacci oscillators through the energy spectrum in the Einstein solid and thus 
expanded to the Debye model, where our results show that the (q { , r/ 2 )-dcformcd Debye temper¬ 
ature, specific heat, thermal and electrical conductivity of ‘deformed’ chemical elements can as¬ 
sume similar values of other pure elements. The results obtained in our study show that by insert¬ 
ing two deformation parameters qi and q 2 , rather than of a parameter q, increases the adjustment 
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range, i.e., we can have different combinations of values as present in Fig. ©. The existence of 
more degrees of freedom as in the present case of two deformation parameters, q\ and q 2 , can be 
well associated with different types of deformations related to two distinct phenomena of disor¬ 
ders or impurities such as, for instance, one due to pressure generating disorders and other due to 
doping, respectively. 
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